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I.  INTRODUCTION 


Traditionally,  physical  optics  and  perturbation  theories  have  been  used  to  derive 
the  like  and  cross  polarized  scattering  cross  sections  for  composite  random  rough 
surfaces.  To  this  end  two-scale  models  have  been  adopted  and  the  rough  surfaces 
are  regarded  as  small  scale  surface  perturbations  that  are  superimposed  on  large 
scale,  filtered  surfaces.  Using  this  model  the  scattering  cross  sections  are  expressed 
as  sums  of  two  cross  sections.  In  general,  however,  the  restrictions  on  both  the  large 
and  small  scale  surfaces  cannot  be  satisfied  simultaneously  using  the  perturbed- 
physical  optics  approaches.  Furthermore  the  results  based  on  the  two-scale  model, 
depend  on  the  wavenumber  where  spectral  splitting  is  assumed  to  occur  between  the 
large  and  small  scale  surfaces.  The  perturbed  physical  optics  approach  (based  on  a 
two  scale  model)  cannot  be  used  to  correctly  predict  the  backscatter  cross  polarized 
cross  sections  since  the  cross  polarized  contribution  from  the  large  scale  surface  is 
assumed  to  be  zero.  Thus,  the  two-scale  model  cannot  be  effectively  applied  to 
polarimetric  radar  studies. 

More  recently,  the  full  wave  solution,  based  on  a  rigorous  mathematical  ap¬ 
proach,  has  been  used  to  determine  the  scattering  cross  sections  for  composite 
random  rough  surfaces  of  finite  conductivity.  Through  the  use  of  rough  surface 
height  conditional  joint  characteristic  functions  (that  account  for  the  rough  surface 
lieight/slope  correlation);  it  is  shown  that  the  same  full  wave  solutions  for  the  dif¬ 
fuse  scattered  fields  are  in  good  agreement  with  the  perturbation  solutions  when 
the  rms  height  and  slopes  are  of  the  same  order  of  smallness  and  with  the  physi¬ 
cal/geometrical  optics  solutions  at  high  frequencies.  The  full  wave  solutions  account 
for  Bragg  scattering  and  specular  point  scattering  in  a  unified  self-consistent  man¬ 
ner.  Using  the  unified  full  wave  approach,  it  is  possible  to  analyze  more  realistic 
models  of  propagation  paths  over  rough  terrain. 
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The  full  wave  approach  can  be  used  to  evaluate  the  bistatic  scattering  cross 
sections  for  all  incident  and  scattered  wave  directions.  These  results  have  been  val¬ 
idated  by  comparisons  with  insitu  measurements  and  laboratory  experiments  and 
numerical  solutions  based  on  Monte  Carlo  simulations  (Bahar  and  Lee  1994).  Data 
on  scattering  from  rough  surfaces  indicate  that  the  full  wave  approach  correctly  pre¬ 
dicts  the  like  to  cross  polarized  scattering  cross  section  ratio  for  rough  seas  (Daley  et 
al  1970,  Bahar  and  Fitzwater  1985).  Furthermore,  the  full  wave  approach  correctly 
predicts  backscatter  enhancement  observed  in  controlled  laboratory  experiments 
using  fabricated  surfaces  with  different  roughness  scales  (Bahar  and  El-Shenawee 
1994).  It  has  also  been  used  to  correctly  interpret  the  Lunar  surface  data  (Bahar 
and  Haugland  1994). 

In  the  work  based  on  the  full  wave  approach,  the  medium  below  the  rough 
interface  is  characterized  by  a  complex  permittivity  to  account  for  its  finite  con¬ 
ductivity.  This  has  an  important  impact  upon  both  single  scattering  and  multiple 
(double)  scattering  by  rough  surfaces.  The  results  for  the  radar  cross  sections  are 
consistent  with  measurements  conducted  in  a  controlled  laboratory  environment. 
The  full  wave  approach  has  also  been  applied  to  rough  surfaces  with  non  Gaussian 
surface  height  probability  density  functions  as  well  as  to  irregular  stratified  media. 
Thus  reflection  and  transmission  across  thin  films  with  rough  interfaces  can  be  in¬ 
vestigated.  The  full  wave  results  are  fully  polarimetric.  These  polarimetric  full 
wrave  investigations  wall  enhance  the  potential  for  remote  sensing  since  they  include 
magnitude  and  relative  phase  data  contained  in  the  sixteen  Mueller  (Stokes)  matrix 
elements. 

II.  MOTIVATION  AND  RESEARCH  OBJECTIVES 

Physical  optics  and  perturbation  theories  have  been  used  to  derive  the  like  and  cross 
polarized  scattering  cross  sections  for  composite  random  rough  surfaces  (Beckmann 
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and  Spizzichino,  1963;  Beckman.  1968;  Rice.  1951).  To  this  end  two-scale  mod¬ 
els  have  been  adopted  and  the  rough  surfaces  are  regarded  as  small  scale  surface 
perturbations  that  are  superimposed  on  large  scale,  filtered  surfaces  (Wright,  1968: 
Valenzuela,  1968;  Barrick  and  Peake,  1968). 

Thus  the  scattering  cross  sections  are  expressed  as  sums  of  two  cross  sections. 
The  first  accounts  for  specular  point  scattering.  It  is  given  by  the  physical  optics 
cross  section  for  the  filtered  surface  consisting  of  the  large  scale  spectral  components. 
The  second  accounts  for  Bragg  scattering.  It  is  given  by  the  cross  section  for  the 
surface  consisting  of  the  small  scale  spectral  components  that  ride  on  the  filtered 
large  scale  surface. 

On  applying  the  perturbed-physical  optics  approaches  it  is  necessary  to  specify 
the  wavenumber  kd  where  spectral  splitting  is  assumed  to  occur  between  the  large 
and  small  scale  spectral  components  of  the  rough  surface.  In  general  the  restrictions 
on  both  the  large  and  small  scale  surfaces  cannot  be  satisfied  simultaneously  and 
using  the  perturbed-physical  optics  approaches  the  evaluation  of  the  scattering  cross 
sections  critically  depends  on  the  specification  of  kd  (Brown.  1978). 

More  recently  the  full  wave  solution  based  on  a  rigorous  mathematical  ap¬ 
proach  has  been  used  to  determine  the  scattering  cross  sections  for  composite  ran¬ 
dom  rough  surfaces  of  finite  conductivity  (Bahar.  1981b:  Bahar  and  Barrick.  1983). 
The  full  wave  solutions  account  for  Bragg  scattering  and  for  specular  point  scat¬ 
tering  in  a  unified  self-consistent  manner.  However,  when  such  a  decomposition  is 
implemented,  the  full  wave  solutions  for  the  scattering  cross  sections  are  expressed 
in  terms  of  a  weighted  sum  of  two  cross  sections  (Bahar.  1981b;  Bahar  and  Barrick, 
1983).  Furthermore,  in  an  attempt  to  draw  more  definite  conclusions  regarding  the 
choice  of  kd,  it  was  varied  over  a  wide  range  of  values  (Bahar  et  al.,  1983a).  It  was 
shown  that  while,  as  expected,  the  like  polarized  cross  sections  associated  with  the 
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large  and  the  small  scale  surfaces  critically  depend  upon  the  choice  of  k^,  the  full 
wave  weighted  sums  of  the  like  polarized  cross  sections  remain  relatively  insensi¬ 
tive  to  variations  in  k^  .  The  physical  optics  approximation  for  the  cross  polarized 
backscatter  cross  section  is  zero.  As  a  result,  the  cross  polarized  backscatter  cross 
section  for  the  large-scale  surface  is  set  equal  to  zero  when  the  two-scale  model  is 
used.  However,  for  backscatter,  only  the  specular  points  on  the  rough  surface  do 
not  depolarize  the  incident  wave. 

It  is  shown  that  while  the  full  wave  solutions  for  the  like  polarized  scattering 
cross  sections  based  on  the  two-scale  model  are  in  reasonable  agreement  (within  3 
db)  with  the  unified  full  wave  solutions,  the  two  solutions  for  the  cross  polarized 
cross  sections  differ  very  significantly  for  near  normal  incidence  (about  15  db). 

In  a  recent  paper  (Bahar,  1991)  the  full  wave  analysis  has  been  provided  for 
rough  surface  diffuse,  incoherent  radar  cross  sections  with  height-slope  correlation 
included.  Since  the  full  wave  rough  surface  element  scattering  coefficients  are  slope 
dependent,  the  expressions  for  the  scattering  cross  sections  involve  the  conditional 
joint  characteristic  functions  for  the  rough  surface  height.  It  is  shown  (Bahar  1991) 
that  when  the  root  mean  square  height  is  small  compared  to  wavelength  and  the 
slope  is  of  the  same  order  of  smallness,  the  full  wave  solutions  reduce  to  Rice:s 
(1951)  small  perturbation  solution.  However,  the  same  full  wave  expression  for  the 
incoherent  scatter  cross  sections  reduces  to  the  physical/geometrical  optics  solution 
in  the  high  frequency  limit  (Beckmann  and  Spizzichino,  1963).  This  is  because,  in 
the  high  frequency  limit,  the  Foureir  transform  of  the  conditional  joint  characteristic 
function  is  proportional  to  the  Dirac  delta  functions  for  the  slopes  at  the  stationary 
phase  (specular)  points  on  the  surface. 

Using  the  unified  full  wave  approach  it  is,  therefore,  possible  to  analyze  more 
realistic  models  of  propagation  paths.  It  is  not  necessary  to  resort  to  the  artificial 
decomposition  of  the  irregular  terrain  into  large  and  small  scale  surfaces. 


7 


The  full  wave  approach  can  be  used  to  evaluate  the  bistatic  scattering  cross 
sections  for  all  incident  and  scattered  wave  directions.  These  results  have  been  vali¬ 
dated  by  comparisons  with  insitu  measurements  and  laboratory  experiments.  Work 
on  scattering  from  rough  surfaces  indicates  that  the  full  wave  approach  correctly 
predicts  the  like  to  cross  polarized  scattering  cross  section  ratio  for  rough  seas  (Da¬ 
ley  et  al..  1970.  Bahar  and  Fitzwater,  1985).  Furthermore,  the  full  wave  approach 
correctly  predicts  backscatter  enhancement  observed  in  controlled  experiments  using 
fabricated  surfaces  with  different  roughness  scales  (Bahar  and  El-Shenawee.  1991a. 
1991b.  1994). 

In  this  work  the  medium  below  the  rough  interface  is  characterized  by  a  com¬ 
plex  permittivity  that  accounts  for  its  finite  conductivity.  This  has  an  important 
impact  upon  both  single  scattering  and  multiple  scattering  by  rough  surfaces.  Work 
on  the  computation  of  the  full  wave  single  and  double  scatter  fields  for  deterministic 
and  random  rough  surfaces  has  been  conducted  (Bahar  and  El-Shenawee.  1991a. 
1991b.  1993).  In  this  work  the  like  polarized  field  scattered  in  the  plane  of  incidence 
is  evaluated  for  both  vertically  and  horizontally  polarized  excitations.  The  effects  of 
varying  the  following  parameters  are  investigated:  angle  of  incidence,  mean  square 
slope,  mean  square  height,  correlation  length  and  complex  permittivity  of  the  scat¬ 
tering  medium.  The  full  wave  results  for  the  enhanced  backscatter  cross  sections 
(due  to  double  scatter)  are  shown  to  be  in  good  agreement  with  experimental  re¬ 
sults  (Knots  et  al.  1993.  Bahar  and  El-Shenawee  1994).  It  is  proposed  to  expand  the 
investigation  of  multiple  scattering  from  two  dimensionally  rough  surfaces.  Thus, 
both  the  like  and  cross  polarized  scattered  field  will  be  considered. 

In  order  to  determine  the  expected  values  of  the  scattering  cross  sections  for 
random  two  dimensional  rough  surfaces,  it  is  necessary  to  characterize  the  surface  by 
conditional  joint  probability  density  functions  for  the  surface  heights  and  slopes  at 
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two  separate  points.  Such  a  probability  density  function  accounts  for  the  correlation 
between  the  slopes  and  heights.  However,  in  order  to  simplify  the  analysis,  the 
correlation  between  the  surface  heights  and  slopes  are  usually  ignored.  Using  the 
full  wave  approach,  the  correlation  between  heights  and  slopes  at  two  distinct  points 
on  the  rough  surface  are  accounted  for  (Bahar,  1991,  Bahar  and  Lee,  1994).  Both 
one  dimensional  and  two  dimensional  rough  surfaces  (characterized  by  their  surface 
height  conditional  joint  probability  density  functions)  are  considered.  This  work 
sheds  light  on  the  source  of  the  discrepancies  between  the  physical  optics  results 
and  the  small  perturbation  results.  Both  like  and  cross  polarized  bistatic  scatter 
cross  sections  are  considered. 

The  first  and  second  order  small  perturbation  solution  (Rice  1951)  as  well  as 
the  physical/geometric  optics  solution  (Beckmann  1968,  Brown  1978)  for  the  cross 
polarized  radar  cross  sections  in  the  plane  of  incidence  (including  backscatter)  are 
zero.  However,  for  two  dimensional  rough  surfaces  the  cross  polarized  cross  section 
is  not  zero,  since  the  normal  to  the  surface  is  not  in  the  plane  of  incidence.  The 
hybrid  solution  based  on  the  two  scale  model  does  account  for  depolarization  in  the 
plane  of  incidence  (including  backscatter)  due  to  the  small  scale  surface  that  rides 
on  the  large  scale  surface  (Valenzuela  1968,  Brown  1978).  However,  using  the  full 
wave  approach,  it  is  shown  that  both  the  large  scale  and  small  scale  two  dimensional 
surfaces  depolarize  the  incident  waves  (Bahar  and  Fitzwater  1984).  This  is  consistent 
with  observations  even  when  the  slopes  are  not  very  large  and  multiple  scatter  is  not 
significant  (Daley  et.  al.  1970).  Since  cross  polarization  is  included  in  the  evaluation 
of  the  Mueller  matrix  elements  for  backscatter,  this  could  have  a  significant  impact 
on  remote  sensing  applications. 

For  the  small  perturbation  solution  (Rice  1951)  and  for  the  physical  optics 
solution,  based  on  the  Kirchhoff  approximations  for  the  surface  fields  (Backmann  and 
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Spizziahino  1963)  the  surface  element  scattering  coefficients  are  independent  of  slope 
while  for  the  full  wave  solutions,  which  are  invariant  to  coordinate  transformations, 
the  surface  element  scattering  coefficients  are  slope  dependent.  For  this  reason  the 
full  wave  solutions  involve  conditional  joint  characteristic  functions  for  the  surface 
heights,  rather  than  the  slope  independent  joint  characteristic  function  (Bahar  1991). 
It  is  this  unique  feature  of  the  full  wave  solution  that  permits  the  use  of  the  same 
expression  to  account  for  Bragg  scatter  from  the  small  scale  spectral  content  of 
the  rough  surfaces  as  well  as  to  account  for  the  specular  contributions  from  the 
neighborhoods  of  the  stationary  phase  points  on  the  rough  surface. 

In  a  recent  work  (Bahar  and  Lee  1994),  a  very  comprehensive  comparison  has 
been  conducted  between  the  full  wave  results  and  the  corresponding  low  frequency 
perturbation  results  and  the  high  frequency  physical/geometrical  optics  results  as 
well  as  to  experimental  data  and  numerical  data  based  on  Monte  Carlo  simulations 
of  rough  surfaces.  It  is  shown  that  when  the  rms  height  and  slope  are  of  the  same 
order  of  smallness,  the  full  wave  results  are  in  excellent  agreement  with  the  small 
perturbation  results  while  at  high  frequencies  they  are  in  very  good  agreement  with 
the  physical/geometrical  optics  results.  There  is  also  good  agreement  with  the 
experimental  and  numerical  results  for  both  polarizations.  This  work  has  also  been 
extended  to  scattering  upon  transmission  across  rough  surfaces. 

Another  important  phase  of  these  investigations  is  related  to  the  effects  of 
shadowing  (Bahar.  1991).  For  example,  a  point  on  the  surface  is  on  the  shadow 
boundary  when  the  scalar  product  of  the  vector  normal  to  the  surface  n  and  the 
incident  wave  vector  in  the  direction  of  propagation  kl  changes  sign.  This  is  referred 
to  as  self  shadow.  A  point  on  the  surface  may  also  be  shadowed  by  a  topographical 
feature  (hill)  not  immediately  adjacent  to  it.  Shadowing  in  this  case  is  character¬ 
ized  by  the  conditional  probability  P2  (Sancer,  1969)  that  a  point  on  the  surface  is 
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both  visible  at  the  receiver  and  illuminated  by  the  source  given  the  rough  surface 
height  and  slope  characteristics.  When  the  full  wave  approach  is  used,  self  shad¬ 
owing  (which  is  accounted  for  intrinsically)  is  shown  to  have  a  significant  effect  on 
scattering.  Using  the  physical  optics  solution  or  the  small  perturbation  solution, 
the  self  shadowing  function  is  ignored.  This  work  on  shadowing  (which  uses  the  full 
wave  approach)  also  has  a  major  impact  on  multiple  (double)  scatter  (Bahar  and 
El-Shenawee  1991a.  b,  1994). 

The  full  wave  approach  has  been  used  to  determine  the  bistatic  reflection  and 
transmission  scattering  cross  sections  of  rough  surfaces  for  all  polarizations,  frequen¬ 
cies  and  incident  and  scatter  wave  directions.  Full  wave  expressions  for  vertically 
and  horizontally  polarized  field  intensities  that  are  diffusely  scattered  across  rough 
interfaces  are  needed  in  order  to  evaluate  the  backscattered  intensities  when  scat¬ 
tering  occurs  at  more  than  one  interface.  The  Principal  Investigator  has  conducted 
extensive  work  in  this  field  (Bahar,  1988;  Bahar  and  Fitzwater,  1989).  This  work 
on  scattering  upon  transmission  across  a  rough  surface  has  been  extended  to  two 
dimensionally  rough  surfaces  and  both  the  like  and  cross  polarized  scattering  cross 
sections  are  accounted  for  in  the  analysis. 

To  enhance  the  potential  applications  of  the  full  wave  approach  to  remote  sens¬ 
ing  problems,  computer  codes  have  been  developed  to  compute  the  sixteen  elements 
of  the  backscatter  Mueller  matrix  (Bahar  and  Kubik.  1989).  This  work  has  been 
used  to  determine  the  tilt  modulation  of  the  scattering  cross  sections  (Bahar  and 
Kubik,  1993a,  1993b).  We  propose  to  generalize  this  phase  of  our  work  in  order  to 
apply  the  full  wave  approach  to  broad  range  of  polarimetric  radars  (Haugland  et  al.. 
1992). 

III.  OVERVIEW  OF  FULL  WAVE  APPROACH 

Rigorous  closed  form  solutions  for  the  reflection  and  transmission  of  electromagnetic 
waves  have  been  derived  for  multilayered  dielectric  structures  of  uniform  thickness 


11 


(Wait,  1962).  However,  in  a  large  variety  of  pertinent  radio  wave  propagation  prob¬ 
lems  the  thicknesses  of  the  layers  are  nonuniform  and  the  height  of  the  interface 
between  two  adjacent  dielectric  layers  is  a  random  variable.  In  these  cases  the  in¬ 
cident  waves  are  depolarized  and  scattered  into  both  propagating  and  evanescent 
waves.  Furthermore,  an  incident  plane  wave  may  be  coupled  into  guided  surface 
waves  and  lateral  waves  of  the  structure. 

■  Often  the  problem  that  is  actually  solved  is  a  highly  idealized  version  of  the 
original  problem  and  concepts  such  as  ”  effective  dielectric  coefficient”  and  ”  effective 
surface  impedance”  are  introduced  in  order  to  make  the  solution  of  the  original 
problem  more  tractable.  However,  the  validity  of  such  approximations  is  limited 
and  they  do  not  necessarily  satisfy  reciprocity. 

Using  a  full  wave  approach,  it  is  possible  to  analyze  more  realistic  models 
of  the  original  physical  structure  without  introducing  simplifying  approximations 
that  cannot  be  justified  apriori  (Bahar  1973a. b,  1974).  The  full  wave  approach 
is  based  on  the  expansion  of  the  electromagnetic  fields  into  a  complete  spectrum 
of  vertically  and  horizontally  polarized  (propagating  and  evanescent)  waves  that 
include  the  radiation  term,  the  lateral  wave  and  the  guided  (surface)  waves  of  the 
structure.  Exact  boundary  conditions  are  imposed  at  each  irregular  interface  of 
the  structure  and  Maxwell’s  equations  are  converted  into  a  rigorous  set  of  coupled 
telegraphists’  equations.  To  this  end.  precise  mathematical  procedures  are  followed 
and  term  by  term  differentiation  of  nonuniform  convergent  series  is  avoided.  The 
coupled  telegraphists’  equations  provide  solutions  for  the  forward  and  backward 
traveling  complex  wave  amplitudes  that  are  related  to  the  electric  and  magnetic 
field  transforms.  Thus,  both  near  field  and  far  field  solutions  can  be  evaluated. 
The  far  fields  scattered  by  the  irregular  layered  structure  can  be  evaluated  in  terms 
of  the  field  transforms  by  using  steepest  descent  techniques.  This  work  has  been 
significantly  refined  in  recent  years  (Bahar  1991,  Bahar  and  Lee  1994). 
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Maxwell's  equations  are  converted  into  the  following  generalized  telegraphists’ 
equations  for  the  forward  and  backward  wave  amplitudes  ap  and  bp  respectively 
(Bahar  1973a.  b). 

-^r-i™P  =  'EY:j{S^a?  +  SBp^)dw'-A’’,  (1) 

Q  v' 

and 

~  +  iubp  =  E  E  /  (Spo«Q  +  s%b°)  dw'  -  Bf.  (2) 

ax  Q  v'  J 

Explicit  closed  form  expressions  for  the  reflection  and  the  transmission  scattering 
coefficients  Spg  have  been  derived  (Bahar  1973b).  Excitations  of  vertically  and 
horizontally  polarized  waves  are  considered.  The  terms  Ap  and  Bp  appearing  in  (1) 
and  (2)  account  for  the  electric  and  magnetic  sources  J,p  and  M,pm. 

The  first-order  iterative  solutions  for  the  wave  amplitudes  ap  and  bp  are  ob¬ 
tained  by  neglecting  the  transmission  and  reflection  scattering  coefficients  in  (1)  and 
(2).  These  first-order  solutions  for  the  primary  fields  are  substituted  on  the  right 
side  of  (1)  and  (2).  and  the  resulting  equations  are  solved  to  obtain  the  second  order 
iterative  solution  for  the  wave  amplitudes.  These  second-order  iterative  solutions 
are  used  in  the  complete  expansions  for  the  electromagnetic  fields  to  obtain  the 
desired  iterative  solutions  for  the  scattered  radiation  fields  through  the  use  of  the 
steepest  descent  method.  Thus  the  first-order  solutions  to  (1)  and  (2)  are  the  pri¬ 
mary  vertically  and  horizontally  polarized  fields  excited  by  the  vertical  electric  and 
magnetic  dipoles  respectively.  The  second-order  iterative  solutions  which  account 
for  depolarization  and  single  scatter  in  arbitrary  directions  are  suitable  even  when 
the  rough  surface  slopes  are  not  small  (Bahar,  1981a).  Since  the  full  wave  expres¬ 
sions  for  the  fields  are  valid  for  all  observation  points,  they  can  also  be  used  to 
determine  the  fields  near  the  surface.  Thus  the  single  scattered  field  that  unpinges 
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upon  the  rough  surface  can  be  used  to  determine  the  double  scattered  field  from 
rough  surfaces  (Bahar  and  El-Shenawee  1994). 

The  principal  results  are  summarized  here  since  some  of  the  recent  references 
are  still  in  press.  The  expression  for  the  vertically  and  horizontally  polarized  diffuse 
single  scattered  fields  from  the  rough  surface  y  =  h{xs)  are  given  by 

EPF  =  G00  J  SFP  [exp(C  •  rs )  -  exp(i^  •  ft)}  dxsEp'  =  EPj  -  EPJ  (3) 

in  which  Spp  is  the  surface  element  scattering  coefficient.  rs  is  the  position  vector  to 
a  point  on  the  rough  surface  and  ft  is  its  projection  on  the  mean  surface  y  =  0.  The 
vector  v  =  k£  —  kl0  where  k£  and  kl0  are  the  scattered  and  incident  wave  vectors,  k0  is 
the  free  space  wave  number  and  G00  is  a  constant.  The  zero  order  (quasi  specular) 
scattered  field  E^  is  not  included  in  the  expression  for  the  diffuse  scattered  field 
(Rice  1951).  Using  (3)  the  corresponding  expression  for  the  ensemble  average  of  the 
scatter  cross  sections  per  unit  area  are  given  by 


<  <j  >= 


Sp®{nf.nl)  Q(n*  ,nl) 


(4) 


For  homogeneous  isotropic  rough  surfaces 

2  id  roo  r  - 

Q(n^ ,  nl)  =  —  /  [X2{vy,-v) -\X(vy)\  ]j0{vxzfd)dfd  (5) 

vy  Jo 

where  =  fti—ft2  X  and  X2  are  the  characteristic  and  joint  characteristic  functions. 
For  the  composite  rough  surfaces  with  large  slopes,  the  scatter  cross  sections  are 
obtained  by  viewing  the  surface  as  a  superposition  of  arbitrarily  oriented  pixels. 
(Bahar  et  al  1983b).  For  each  arbitrarily  oriented  pixel,  the  scatter  cross  section  is 
given  by 

(<rp«)p  =  |l3p«j!Q!.(«/,ni)  (6) 
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where  DPQ  is  the  element  of  the  scattering  matrix  D  given  by 


D  =  TfSpr 


(7) 


in  which  Sp®,  the  elements  of  the  2x2  matrix  Sp  is  obtained  from  Sp®  on  replacing 
the  angles  of  incidence  and  scatter  in  the  reference  coordinate  system  by  the  inci¬ 
dent  and  scatter  angles  in  the  (local)  pixel  coordinate  system  and  the  2x2  matrices 
Tl  and  T 7  relate  the  vertically  and  horizontally  polarized  waves  of  the  reference 
coordinate  system  (primed)  to  the  vertically  and  horizontally  polarized  waves  of  the 
pixel  (unprimed)  coordinate  system.  Furthermore. 


9Z'-* 

=  0 


[X2(Vy,  ~Vy)  ~  \X(Vy)f  JQ  (V  XZT  <f)  f  d 


(8) 


K  •  avYvl Jo 

The  scatter  cross  section  for  the  composite  rough  surface  is  obtained  by  summing 
the  radar  cross  sections  of  the  individual  pixels.  Thus  if  the  slopes  of  each  pixel  is 
also  random  (with  probability  density  function  p(hx,  hz).  the  scatter  cross  section 
(per  unit  area)  for  the  composite  rough  surface  (aPQ)c  is  given  by  the  statistical 
average  of  the  pixel  scatter  cross  section  per  unit  area 


(aPQ  =  (\DPQ\2QP(nf:ni) 

The  surface  height  autocorrelation  function  for  the  pixel  is 

(hsu  hS2 )  =  (h2s)Rs(rd)  =  2tt  ]-^J0(krd)kdk 

Jkd  4 

and  the  total  mean  square  slope  of  the  pixel  is 

4  =  2 

Jo  4 


(9) 


(10) 


(11) 


where  kd  =  2i r/LP  (LP  is  the  lateral  dimension  of  the  pixel)  kc  is  the  cut  off  wave 
number  (Brown,  1978)  and  W(k)  is  the  surface  height  spectral  density  function. 
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It  is  readily  shown  that  when  the  surface  heights  and  slopes  are  of  the  same 
order  of  smallness  the  full  wave  solution  (4)  reduces  to  the  small  perturbation  solu¬ 
tion. 

A  third  order  iterative  solution  to  the  generalized  telegraphists  equations  (1). 
(2)  has  also  been  carried  out  in  order  to  obtain  the  full  wave  results  for  the  multiple 
(double)  scattered  field  (Bahar  and  El-Shenawee  1993,  1994). 

IV.  ENHANCED  BACKSCATTER 

Enhanced  backscatter  from  rough  surfaces  has  been  observed  in  numerous  carefully 
conducted  experiments.  While  the  original  measurements  of  backscatter  enhance¬ 
ment  were  conducted  using  two  dimensionally  rough  surfaces,  no  numerical  or  an¬ 
alytical  solutions  for  the  scattered  electromagnetic  fields  were  available.  Thus  it 
was  not  possible  to  conduct  a  parametric  study  of  the  observed  phenomena  and  to 
physically  interpret  them.  More  recently,  experimental,  numerical  and  analytical 
investigations  have  been  carried  out  on  surfaces  that  are  essentially  rough  in  only 
one  dimension.  The  plane  of  incidence  was  restricted  to  the  plane  in  which  the  local 
normal  to  the  rough  surfaces  lies.  Thus  these  surfaces  do  not  depolarize  the  incident 
waves  and  the  problem  is  essentially  scalarized.  While  these  one  dimensional  rough 
surfaces  do  not  have  the  same  range  of  physical  applications  as  the  two  dimen¬ 
sional  rough  surfaces,  the  numerical  and  analytical  solutions  to  the  electromagnetic 
scattering  problem  for  the  one  dimensional  rough  surfaces  are  easier  to  interpret. 

This  intermediate  step  provides  the  necessary  impetus  to  return  to  the  hitherto 
intractable  (both  numerically  and  analytically)  problem  of  electromagnetic  scatter¬ 
ing  from  two  dimensionally  rough  surfaces.  To  date  no  comprehensive  computer 
code  can  be  used  to  solve  the  problem  of  scattering  from  two  dimensionally  rough 
surfaces  numerically  (moment  methods)  in  a  tractable  manner.  Computers  are  used 
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to  conduct  a  complete  parametric  study  of  the  scattering  phenomena.  This  para¬ 
metric  study  involves  the  determination  of  the  effects  of  changing  the  statistical 
characterizations  of  the  rough  surface  upon  the  enhanced  backscatter  phenomena. 
The  random  rough  surface  is  characterized  by  the  joint  surface  height/slope  prob¬ 
ability  density  function  at  two  points  on  the  surface  (including  the  mean  square 
height  and  slope,  the  surface  height  autocorrelation  function,  and  the  correlation 
length).  The  effects  of  changing  the  incident  angle,  the  polarization  of  the  incident 
and  scattered  waves,  and  the  medium  parameters  are  studied.  The  research  was 
initially  restricted  to  the  derivation  of  high  frequency  physical  optics  approxima¬ 
tions  since  it  sheds  light  on  the  analytical  and  numerical  solutions  to  the  problem. 
The  precise  level  of  the  polarimetrie  backscattered  intensities  (the  Stokes  vectors 
related  by  the  Mueller  matrix),  as  well  as  the  range  of  angles  about  the  backscatter 
direction  for  which  the  observed  enhanced  backscattered  is  significant  is  studied  in 
detail. 

V.  OUTLINE  OF  APPROACH  FOR  DOUBLE  SCATTER  CROSS  SEC¬ 
TIONS  AND  SUMMARY  OF  MOST  IMPORTANT  RESULTS 

The  expressions  for  the  vertically  and  horizontally  polarized  diffuse  single  scat¬ 
tered  radiation  fields  from  the  two-dimensionally  rough  surface  ys  =  h(xs.  zs){—L  < 
xs,  zs  <  L)  are  given  in  matrix  form  as  follows: 

G{(f)  =  ( A)2  /  /  S{i,  i'JexpH*'  -  f) 

exp  (iv1  ■  rs  —  exp  (iv'  ■  ft) 

-  ^ 

in  which  G{  and  Gl  are  2  x  1  matrices  whoe  elements  are  the  vertically  and 
horizontally  polarized  components  of  the  scattered  and  incident  fields  Evi  and 
EPl(P  =  V,  H)  respectively.  The  elements  of  the  2x2  scattering  matrix  S(k'.kl) 
are  the  like-and  cross-polarized  scattering  coefficients  SFQ(P,Q  =  V,  H).  These 


dXsdzJkkdk'^ 

S  S  IJ 
VOX 


G^o) 


(12) 
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coefficients  are  functions  of  the  incident  and  scatter  wave  vectors  kl0  —  kQnl  and 
k'0  —  k0n'  respectively  and  the  electromagnetic  parameters  of  the  media  on  both 
sides  of  the  rough  interface  (Bahar  and  Lee  1994).  The  position  vectors  to  the  ob¬ 
servation  point  and  to  a  point  on  the  rough  surface  are  f  and  fs  respectively  and  ft 
is  the  projection  of  fs  on  the  mean  plane  y  =  0.  The  vector  v  is  given 


=  k' 


K  =  vx 


dy  T  VyCLy  "f"  VZ 


(13) 


The  integration  in  (3)  is  over  the  rough  surface  variables  xs  and  zs  and  the  wave 
vector  variables  k'oz  and  k'  .  In  (3)  exp (iu>t)  time  excitations  are  assumed.  Note 
that  the  expression  for  G{  remains  finite  as  v'y  — >0  (grazing  incident  and  scatter 
angles) . 

For  simplicity  we  initially  consider  one-dimensional  rough  surfaces  y  =  h(xs). 
The  incident  and  wave  vectors  lie  in  the  x,  y  plane  and  there  is  no  depolarization 
Sp®(k,  kl )  — >  0(P  ^  Q).  In  this  case  the  vertically  and  horizontally  polarized  waves 
are  uncoupled  and  the  scattering  problem  is  scalarized.  Thus  equation  (12)  can  be 
expressed  as  follows: 


GfsP(r)  =  ^[  I  Sp{k'P)exp(-ik'o:kl0) 

Z7T  Jn'=— oo  Jxs=—L 

exp (iv1  ■  fs)  —  exp(zu'  •  f,)  —  exp (iv1  ■  ft) 


dk' 


U{fs)dx'-^Gip(o).  P  —  V.  H(14) 

rC 


in  which  U(fs )  is  a  shadow  function  (Sancer  1969).  Equation  (14)  is  used  to  obtain 
the  differential  wave  vector  contribution  to  the  diffuse  single  scatter  field  that  im¬ 
pinges  on  the  surface  at  fs2(rs  — *  fsl)  and  (f  — »  fs2).  This  expression  is  used  in  (14) 
instead  of  the  incident  plane  wave  GlP  and  the  expression  for  the  double  scattered 
field  is  obtained  (Bahar  and  El-Shenawee  1944).  Use  is  made  of  the  steepest  de¬ 
scent  method  to  obtain  the  far  field  expression  for  the  double  scattered  field  in  the 
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direction  of  the  wave  vector  k{  =  k0n^ . 


=  j\  J'  Sp(V,n') 

v0T  Jn'=- 1  Jx'=-L 


47 T2  v  k 

exp {-ik0x'sl(nlx  -  n'x)}  exp {ik0x's2(nfx  -  n'x )} 
exp{*&0h(:r(2)(n^  -  nj,)}  -  1 


x 


[exp(-ik0h(x'sl)(ni  -  i 

<j))  -  l] 

ko(~ni  +  n'y) 

x 


dn[ 


ko(n(  -  n'y) 


U{f'sl)U{fs2)dx'sldx's7-^G'p{  0) 

77-„ 


(15) 


The  expression  for  the  double  scattered  field  (14)  is  multiplied  by  its  complex  con¬ 
jugate  to  obtain  the  expression  for  the  double  scatter  power.  After  tedious  manipu¬ 
lation  of  the  expression  for  the  scattered  power  (and  accounting  for  the  height-slope 
correlations  [Bahar  19911)  it  is  shown  that  the  (mean)diffuse  scatter  cross  section 
per  unit  area  can  be  expressed  as  follows  for  the  quasi  parallel  scattering  path 
n'x  =  yfi  ~  ny,nx  =  \Jt  ~  n'y 2  and  n'x  =  -yjl  -  n*,n"  =  -l^Jl  -  n"2  (Bahar  and 
El-Shenawee  1994). 

2  k0Lm 

\tism)rz\n- )r2\n‘ ) 


{a pd. )  — 


27 r2 


:(2 Im)P2(hi)P2(h/)  /  [  {vpdi{n'y.  n"y)(apd2(n'y,  n"y)) 

”  71  ftu  J  Tin  ii 


x [1  P2(  n'y  )]smc{k0Lmndx}exp(-k20(h2)n2dy-  dHay 


dn: 


dy 


(16) 


•J1  -  K)!  \jli  -  K)2 

where  Lm  is  the  width  of  a  typical  depression  on  the  rough  surface.  <  apdq{n' ,  n")  > 
and  <  <jpd2(n'y,n'y )  are  associated  with  single  scatter  cross  sections  and  P2  are  the 
shadow  functions  (Sancer  1969).  The  vectors  na  and  nd  are  defined  as 


n  4-  n" 


nn  = 


-  -  I  -  // 

rid  “  n  —  n 


(17) 


The  major  contribution  to  the  integral  (16)  is  around  nd  =  O  for  the  quasi  parallel 
case  due  to  the  interaction  terms  (between  the  single  scatter  cross  sections  <  a  >pdi 
and  <  apd2  >)  sinc(k0Lmndx)exp(-k20  <  h 2  >  n2dy)  in  the  integral  (16). 

The  expression  for  the  diffuse  double  scatter  cross  section  associated  with  the 
quasi  antiparallel  paths  (n'x  =  ^1  -  ra^.n"  =  -yj  1  -  n"2andn(.  =  -yj  1  -  n"  = 
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yj  1  —  n"2)  can  be  expressed  as  follows: 

{(Tad)  =  P2(nl)P2{nf)  I  j  (aadl{n'y,  n"y)){aad2(n'y)  n"y))x 

£i/T  J  TIq  y  JTlay 

t1  -  P 2 (l^y D] [1  ~  P2(\n'y\)]smc[k0Lm{nfx  +nlx-  2 nax)} 

In  (18)  {(Jadiin'y.n'y))  and  (aad2(n'y.n"y))  are  associated  with  single  scatter  cross 
sections.  In  (16)  and  (17).  the  region  of  integration  is  half  the  diamond  shaped  area 
with  end  points  on  the  nay  axis  at  (—1, 1)  and  on  the  ndy  axis  at  (-2,  2).  The  major 
contributions  to  the  integral  (18)  is  around  na  =  O.  Due  to  the  interaction  terms 
sinc[k0Lm(nl  +  ntx-2nax)}  exp (—kl(h2)(n^  +  ny  -  2nfly)2]  there  is  a  sharp  peak  in  the 
expression  for  (a pd)  for  the  backscatter  direction  hf  =  -nl  (enhanced  backscatter). 

A  comprehensive  parametric  study  has  been  conducted  to  determine  the  de¬ 
pendence  of  the  level  and  the  angular  width  of  the  sharp  enhanced  backscatter  upon 
the  rough  surface  characteristics  (mean  square  height  (h2)  and  slope  (h2x)  and  the 
mean  width  of  a  typical  depression  Lm)  (Bahar  and  El-Shenawee  1994).  A  high  fre¬ 
quency  stationary  phase  (specular  point  approximation  of  the  full  wave  results)  has 
also  been  obtained.  These  significantly  simplified  expressions  shed  light  on  the  phys¬ 
ical  interpretation  and  the  numerical  evaluation  of  the  double  scatter  cross  section. 
However,  it  is  shown  that  these  high  frequency  results  are  practically  independent 
of  polarization  for  highly  conducting  surfaces.  The  angular  width  of  the  enhanced 
backscatter  cross  section  is  also  narrower  when  the  high  frequency  approximations 
are  used. 


We  now  turn  to  multiple  scatter  from  the  more  general  and  interesting  two 
dimensionally  rough  surfaces.  Proceeding  as  in  the  case  of  the  one  dimensional  rough 
surface,  the  expression  for  the  double  scattered  far  field  is  expressed  as  follows: 

GfsP(rs )  =  J  S2'{nf ,  nl)exp(ik0nf  ■  rs2)exp[-ikl  ■  fj 


Sj'(n',  n*)exp(-?T0nz  ■  r'sl)- 


dn'ydn'z 


--dxsXdz[1dx's2dzls2G1{o) 
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The  expression  for  the  diffuse  double  scattered  intensity  is  obtained  on  multiplying 
(19)  by  its  complex  conjugate.  Thus  the  diffuse  double  scatter  cross  sections  (crpd) 
and  (<7ad)  associated  with  the  quasi  parallel  and  quasi  anti-parallel  paths  respectively 
can  be  expressed  as  follows: 

M  =  8(2)°jr-)2  /  /  /  /(^i)(^2)[i-J32(KI)][i-P2(KI)] 

smc(k0Lmndx)smc(k0Lmndz)exp(-k20(h2)n2dy 

dnaydTlazdTldy(lTldz  tool 

(1  -  n'2  -  n,2){{  1  -  n"2  -  n"2)]a 

in  which  Lm  is  the  mean  width  of  a  typical  depression  on  the  surface  and  (apdi)  and 
(<7pd2)  are  associated  with  single  scatter  cross  sections  for  the  quasi  parallel  paths. 
Furthermore 


Wad)  = 


8(2^m)i  / // /(^>(^2)[l-P2(KI)](l-P2(KI)] 

sine [k0Lm{n{  +  n\  -  2nax)]sinc[k0Lm(nfz  +  n\  -  2naz)}exp[-k70(h 2 


K  +  <  -  2nly)\- 


dnay  dnaz  dvidydTidz 


— T  (21) 

(1  -  n'y  -  nj?)[(  1  -  n"2 0  -  n"2)]2 
In  (21)  (cradl)  and  (aad2)  are  associated  with  single  scatter  cross  sections  for 
the  quasi  antiparallel  paths.  In  (20)  and  (21)  the  regions  of  integration  are  half  the 
diamond  shaped  areas  with  end  points  on  the  nay  axis  at  (-1,  1)  and  on  the  ndy  axis 
at  (-2.  2)  and  on  the  naz  axis  at  (-1,  1)  and  on  the  ndz  axis  at  (-2,  2)  (Bahar  and 
El-Shenawee,  1994). 

In  equation  (20)  in  view  of  the  interaction  terms  between  the  single  scatter 
cross  section,  sine (k0Lmndx)  smc(k0Lmndz)exp[-kl(h2)n2dy}.  the  major  contribution 
to  the  integral  is  around  nd  =  O  for  the  quasi  parallel  paths.  Similarly  in  equation 
(21)  in  view  of  the  interaction  terms  between  the  single  scatter  cross  sections  sine 
{k0Lm(nfx  +  nlx-  2 nax)]  sine  [k0Lm(nfz  +  n*  -  2nQZ]  exp  \-k20,  (h2)(nf  +  n*  -  2 nay)2] 
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and  the  shadow  functions  jl  —  Pi(\nfy\)]  [1  —  ln//j/l)]  the  major  contribution  to 

the  integral  is  around  na  =  0  for  the  quasi  antiparallel  paths. 

Using  the  high  frequency  approximations,  the  cross  sections  associated  with 
the  single  scatter  (<rpdl).  (apd 2).  (cradi)  and  (fxad2)  can  be  expressed  in  closed  form 
and  the  resulting  expressions  for  the  double  scatter  cross  section  can  be  expressed 
in  terms  of  four  dimensional  integrals.  Developing  an  efficient  computer  code  for 
equations  (20)  and  (21)  in  the  high  frequency  limit  is  in  progress.  The  results 
sheds  light  on  the  physical  interpretation  and  on  the  numerical  evaluation  of  the 
full  wave  double  scatter  cross  section  in  a  tractable  manner  as  in  the  case  of  the  one 
dimensionally  rough  surface  (Bahar  and  El-Shenawee.  1994).  See  page  21a. 
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X.  CONCLUDING  REMARKS 

The  research  group  at  the  University  of  Nebraska-Lincoln  has  over  the  three  year 
duration  of  the  ARO  Grant,  satisfied  all  the  proposed  research  objectives  (including 
comparisons  with  theoretical,  numerical  and  experimental  results).  The  problems 
of  multiple  (double)  scatter  and  enhanced  backscatter  has  been  analyzed  and  the 
radar  cross  sections  have  been  evaluated  in  a  tractable  manner.  The  effects  of  the 
rough  surface  statical  parameters  on  the  angular  width  and  level  of  the  enhanced 
backscatter  peak  has  been  considered  in  detail. 

Several  new  objectives  need  to  be  considered  for  future  work  in  this  field.  The 
extension  of  this  work  on  double  scatter  from  two  dimensionally  rough  surfaces  need 
to  be  continued  (outline  of  this  work  is  given  in  Section  V).  The  analytical  solutions 
explicitly  relate  the  rough  surface  parameters  (mean  square  slope  and  height,  and 
mean  width  of  typical  depressions  on  the  surface)  to  the  quasi  parallel  and  anti¬ 
parallel  contribution  to  the  double  scatter  cross  sections.  It  is  also  necessary  to 
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exploit  recent  advances  in  polarimetric  radar  techniques  by  measuring  all  sixteen 
Mueller  matrix  elements. 

Applications  to  irregular  stratified  media  should  be  considered.  This  work 
would  also  be  applied  to  Ellipsometry  (quasi  specular  reflection).  Instead  of  em¬ 
ploying  effective  media  parameters  to  model  the  effects  of  rough  surfaces  on  the 
ellipsometric  parameters,  a  more  realistic  physical  model  of  the  rough  surface  ef¬ 
fects  could  be  developed  using  the  full  wave  analysis.  Thus  both  the  coherent  and 
incoherent  components  of  the  Scattering  Cross  Sections  should  be  evaluated. 

A  major  effort  can  now  be  made  at  the  University  of  Nebraska  to  compare 
analytical/numerical  results  with  both  in  situ  measurements  as  well  as  controlled 
laboratory  experiments  using  a  4tt  Optical  Scatterometer  (Reflection  and  Transmis¬ 
sion)  and  Atomic  Force/Scanning  Tunneling  Microscopes.  The  principal  investiga¬ 
tor  will  be  submitting  proposals  to  ARO  to  support  the  proposed  extensions  to  our 
investigations. 
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for  review  (Preprint  enclosed)  (see  Item  #3).  Four  (4)  papers  were 
published  in  journals  (see  Item  #7)  and  reprints  of  these  publications 
have  been  submitted  with  the  report.  The  Ph.D.  dissertation  of.  Ms .  M. 
El-Shenawee  and  the  M.S.  thesis  of  B.  S.  Lee  are  accepted  (see 
Item  #8 ) . 

In  the  paper,  "Full  Wave  Analysis  for  Rough  Surface  Diffuse 
Incoherent  Radar  Cross  Sections  with  Height-Slope  Correlations  Included" 
one  dimensionally  rough  surfaces  are  characterized  by  four  dimensional 
joint  probability  density  functions  for  the  heights  and  slopes  at  two 
points  on  the  surface  at  a  distance  xd  apart.  It  is  shown  that  the 
expression  for  the  full  wave  incoherent  diffuse  scatter  cross  section 
reduces  to  the  low  frequency  small  perturbation  solution  when  the 
surface  heights  and  slopes  are  small  and  to  the  high  frequency  physical 
optics  solution  when  the  major  cocntributions  to  the  scattered  fields 
come  from  the  neighborhoods  of  the  stationary  phase  (specular)  points. 

Using  the  super-computer  facility  at  Cornell  University  (NSF 
sponsored)  the  full  wave  approach  has  also  been  used  to  evaluate 
singly  and  multiply  scattered  electromagnetic  fields  from  rough 
surfaces  (see  Item  #7).  It  is  shown  that  for  highly  reflective 
surfaces,  multiply  scattered  fields  can  make  significant  contributions 
to  enhanced  backscatter  for  near  normal  incidence  if  the  mean  square 
slopes  are  large  (  <h  >-»l)  and  the  correlation  lengths  are  small 

X 

^c/  A^l,  A  is  the  free  space  wavelength. 


The  principal  investigator  has  been  invited  by  Dr.  W.  Flood, 
the  project  monitor,  to  participate  in  a  workshop  on  enhanced 
backscatter  at  the  University  of  Colorado,  at  Boulder,  on  January  10, 
1992  after  the  International  Radio  Science  Meeting.  At  the  Meeting 
researchers  will  present  experimental,  numerical  and  analytical 
studies  of  rough  surface  scattering. 

An  application  has  been  made  to  the  project  monitor  for  matching 
funds  to  upgrade  our  Scanning  Tunnelling  Microscope/Atomic  Force 
Microscope  (STM/AFM)  and  our  optical  polarimetric  scatterometer . 

These  upgrades  should  enhance  our  ability  to  scan  larger  areas  of 
the  rough  surface  with  the  STM/AFM  and  to  enhance  the  performance 
of  the  scatterometer  in  the  backscatter  direction. 
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1.  Papers  Presented  at  Technical  Meetings 

(i)  International  Union  of  Radio  Science  (URSI)  Meeting,  University  of  Colorado,  Boulder, 
CO,  January  7-11,  1992,  "Diffuse  Incoherent  Radar  Group  Cross  Sections  for  Random 
Rough  Surfaces  with  Height-Slope  Correlations ,  Included." 

(ii)  Specialist  Meeting  on  Microwave  Radiometry  and  Remote  Sensing  Applications,  Boulder, 

CO, January  14-16,  1992,  "Remote  Sensing  of  Monomolecular  Oil  Slicks  Characterized  by 
Rough  Surface  Spectra  Obtained  in  Controlled  Experiments,  Unified  Full  Wave  Approach," 
with  R.  D.  Kubik. 

(iii)  1992  International  Geoscience  and  Remote  Sensing  Symposium  (IGARSS'92)  and  International 
Radio  Science  Meeting,  Houston,  TX,  May  26-29,  1992,  "Unified  Full  Wave  Solution  for 
the  Scattering  Cross  Sections  of  the  Rough  Sea  Surface  Coated  with  a  Thin  Film," 

with  R.  D.  Kubik. 

(iv)  Advisory  Group  on  Aerospace  Research  and  Development  (AGARD),  Electromagnetic  Wave 
Propagation  Panel,  50th  Symposium  on  Radio  Location  Techniques,  June  1-5,  1992, 
"Irregular  Media  Effects  on  Radiowave  Signals  Used  in  Navigation  and  Positioning 
Systems." 

(v)  1992  CRDEC  Scientific  Conference  on  Obscuration  and  Aerosol  Research,  Aberdeen,  MD, 

|  June  22-25,  1992,  "Use  of  a  New  Polarimetric  Optical  Bistatic  Scatterometer  to 

|  Measure  the  Transmission  and  Reflection  Mueller  Matrix  for  Arbitrary  Incident  and 

•  Scatter  Direction." 

1 2 .  Papers  Accepted  for  Presentation  at  Technical  Meetings 

1  (i)  1992  IEEE-APS  International  Symposium  and  URSI  Radio  Science  Meeting,  July  20-24, 

1 992 ,  Chicago,  IL,  "Multiple  Scattering  from  Random  Distribution  of  Individual  Rough 
Surface  Scatterers,"  with  M.  El-Shenawee. 

(ii)  URSI  International  Symposium  on  Electromagnetic  Theory,  Sydney,  Australia,  Auqust  17-20, 
1992,  "Transformation  of  Rice's  Small  Perturbation  Results  for  Rough  Surface  Scattering 
into  a  Comprehensive  Single  Scatter  Solution  that  Includes  Physical  and  Geometrical 
Optics." 

3.  Papers  Submitted  to  Journal  Editors  for  Review 

[  0)  "Unified  Full  Wave  Solutions  to  Interpret  Apollo  Lunar  Surface  Data,"  with  M. 

Haugland. 

(ii)  "Full  Wave  Single  and  Multiple  Scattering  from  Rough  Surfaces,"  with  M.  El-Shenawee. 

(iii)  "Use  of  a  New  Polarimetric  Optical  Bistatic  Scatterometer  to  Measure  the  Transmission 
and  Reflection  Mueller  Matrix  for  Arbitrary  Incident  and  Scatter  Directions,"  with 

R.  D.  Kubik  and  D.  R.  Alexander. 

4.  Papers  Accepted  for  Publication 

(i)  "Simulation  of  High  Resolution  Radar  Polarimetric  Images  -  Unified  Full  Wave  Approach," 
with  R.  Kubik,  International  Journal  of  Remote  Sensing. 

(ii)  "Identification  of  Contaminant  Coatings  Over  Rough  Surfaces  Using  Polarized  IR 
Scattering,"  with  S.  M.  Haugland  and  A.  H.  Carrieri,  Applied  Optics. 

(iii)  "Active  Remote  Sensing  of  the  Ionosphere  and  the  Earth  Surface  Irregularities," 
Proceedings  of  the  49th  AGARD-EPP  Symposium  on  Remote  Sensing  of  the  Propagation 
Environment,  Cesme,  Turkey,  1991. 

(iv)  "Remote  Sensing  of  Monomolecular  Oil  Slicks  Characterized  by  Rough  Surface  Spectra 
Obtained  in  Controlled  Experiments:  Unified  Full  Wave  Approach,"  with  R.  D.  Kubik, 
Proceedings  of  the  Specialist  Meeting  on  Microwave  Radiometry  and  Remote  Sensing 
Appl ications,  July  1992. 


(v)  "Multiple  Scattering  from  Random  Distribution  of  Individual  Rough  Surface 
Scatterers,"  with  M.  El-Shenawee,  Proceedings  of  the  1992  Antennas  and 
Propagation  Society  International  Symposium,  July  1992. 

( v i )  "Irregular  Media  Effects  on  Radiowave  Signals  Used  in  Navigation  and 
Positioning  Signals,"  Proceedings  of  the  Advisory  Group  on  Aerospace  Research 
and  Development  (AGARD)  Electromagnetic  Wave  Propagation  Panel,  50th 
Symposium  on  Radio  Location  Techniques,  July  1992. 

(vii)  "Transformation  of  Rice's  Small  Perturbation  Results  for  Rough  Surface 
Scattering  into  a  Comprehensive  Single  Scatter  Solution  that  Includes 
Physical  and  Geometrical  Optics,"  Proceedings  of  the  1992  International  Union 
of  Radio  Science  (URSI)  Symposium  on  Electromagnetic  Theory,  Sydney, 
Australia,  August  1992. 

5.  Papers  Published 

(i)  "Unified  Full  Wave  Solution  for  the  Scattering  Cross  Sections  of  the  Rough 
Sea  Surface  Coated  with  a  Thin  Film  of  Oil,"  with  R.  D.  Kubik,  Proceedings 
of  the  1992  International  Geoscience  and  Remote  Sensing  Symposium  (IGARSS'92) 
and  International  Radio  Science  Meeting,  June  1992. 

6.  Participation  in  Workshop 

Workshop  on  Enhanced  Backscatter  organized  by  W.  Flood,  Army  Research 
Office  (P.0.  Box  12211,  Research  Triangle  Park,  NC)  at  Boulder,  CO, 

January  10-11,  1 992. 


BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


During  this  reporting  period  (January  1,  1992  -  June  30,  1992)  the  Principal  Investigator  par¬ 
ticipated  in  five  (5)  Scientific/Technical  Meetings  where  five  papers  were  presented  (See  Item  #1), 
and  two  (2)  papers  have  been  accepted  for  presentation  at  International  Conferences  during  the 
next  few  months  (See  Item  #2).  Three  (3)  papers  were  submitted  to  Journal  Editors  for  review 
(See  Item  #3).  Seven  (7)  papers  were  accepted  for  publication  in  Technical/Scientific  Journals 
and  Conference  Proceedings  (See  Item  #4).  One  (1)  paper  was  published  (See  Item  #5,  Reprints 
attached).  The  Principal  Investigator  participated  in  the  Workshop  on  Enhanced  Backscatter  in 
Boulder,  Colorado  January  10-11,  1992  organized  by  Dr.  W.  Flood  of  the  Army  Research  Office 
(See  Item  #6). 

The  investigations  on  the  “Full  Wave  Analysis  for  Rough  Surface  Diffuse  Incoherent  Radar 
Cross  Sections  with  Height-Slope  Correlations  Included,”  (See  Semi-Annual  Report  #1)  have  been 
extended  to  vertically  polarized  excitations.  It  is  shown  that  for  random  rough  surfaces  with  very 
small  Rayleigh  roughness  parameters  (3  =  4 k%  <  h2  >  (where  kQ  is  the  free  space  wave  number 
and  <  h2  >  is  the  man  square  height)  and  small  slopes  (hx  «  0.1)  the  full  wave  solutions  for 
the  diffuse  scattered  vertically  polarized  waves  are  in  good  agreement  with  the  small  perturbation 
solution  of  Rice.  As  the  slope  hx  becomes  extremely  small,  the  second  order  small  perturbation 
solution  is  not  valid.  The  full  wave  solution  is  also  compared  with  the  physical  optics  solution. 
The  scattered  vertically  polarized  full  wave  radiation  fields  vanish  at  grazing  angles  in  agreement 
with  experimental  data.  The  corresponding  small  perturbation  and  physical  optics  solutions  do 
not  vanish  at  grazing  angles.  This  work  will  be  submitted  to  a  technical/scientific  journal  for 
publication  during  the  next  reporting  period. 

Progress  has  also  been  made  on  obtaining  solutions  for  multiple  scattering  from  rough  surfaces. 
The  work  based  on  a  model  consisting  of  random  distributions  of  surface  scattering  will  be  presented 
in  July  1992  (See  Item  #2,  (i))  at  the  IEEE-APS  International  Symposium  in  Chicago,  Illinois. 
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Vol .  31,  No.  19,  July  1,  1992. 

(ii)  "Remote  Sensing  of  Monomolecular  Oil  Slicks  Characterized  by  Rough  Surface 
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with  R.  D.  Kubik,  Proceedings  of  the  Specialist  Meeting  on  Microwave 
Radiometry  and  Remote  Sensing  Applications,  June  1992. 
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Propagation  Society  International  Symposium,  July  1992. 
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Brief  Outline  of  Research  Findings 


During  the  reporting  period  (1  July  1992  -  31  December  1992)  the 
principal  investigator  submitted  papers  at  four  (4)  Scientif¬ 
ic/Technical  meetings  (see  Item  #1)  and  ten  (10)  papers  were 
submitted  for  presentation  at  international  conferences  during  1993 
(see  Item  #2) .  Five  (5)  papers  were  submitted  to  Journal  Editors 
for  review  (see  Item  #3).  Three  (3)  papers  were  accepted  for 
publication  in  Technical/ Scientific  Journals  (see  Item  #4)  and  four 
(4)  were  published  in  Technical/Scientific  Journals  and  Conference 
proceedings  (see  Item  #5) .  Reprints  of  Item  #5  are  enclosed.  Also 
enclosed  are  preprints  of  Item  #2  (iii) ,  (v) ,  (ix).  The  topic  in 
these  three  preprints  is  enhanced  backscatter  due  to  multiple 
scatter  from  rough  surfaces. 

Very  significant  progress  has  been  made  in  developing  and  executing 
the  computer  codes  for  the  bistatic  multiple  (double)  scatter  cross 
sections  for  one  dimensionally  random  rough  surfaces.  Both  the 
vertically  and  horizontally  polarized  waves  are  considered.  These 
full  wave  results  are  expressed  in  terms  of  six  dimensional 
integrals  involving  two  spatial  variables,  two  slope  variables,  and 
two  wave  vector  variables.  The  major  contributions  come  from  the 
quasi-parallel  and  quasi  anti-parallel  paths  for  the  field  and  its 
complex  conjugate.  The  contribution  from  quasi  anti-parallel  paths 
which  is  responsible  for  the  enhanced  backscatter  is  also  responsi¬ 
ble  for  the  observed  undulations  in  the  scattered  intensity  on 
either  side  of  the  backscatter  direction.  This  is  due  to  the 
interference  of  the  contributions  from  two  adjacent  regions  of  the 
rough  surface. 

A  high  frequency  approximation  of  the  full  wave  results  reduces  the 
six  dimensional  multiple  scatter  integral  into  a  two-dimensional 
integral  involving  the  wave  vector  variables  only.  These  high 
frequency  approximations  provide  physical  insight  to  the  scattering 
problem.  However,  they  are  practically  independent  of  polariza¬ 
tion.  The  results  obtained  from  the  six  dimensional  integrals  show 
significant  polarization  dependence  of  the  multiply  scattered  field 
even  though  the  single  scattered  field  is  practically  independent 
of  polarization.  The  full  wave  expressions  offer  an  explanation  to 
these  results.  The  three  enclosed  pre-prints  provide  additional 
details  regarding  this  phase  of  our  research. 

Progress  has  also  been  made  on  the  comparison  of  the  numerical 
solutions  to  the  Generalized  Telegraphists'  Equations  with  two 
single  scatter  interactive  solutions.  The  interactive  solutions 
are  based  on  the  use  of  "local"  basis  functions  and  basis  functions 
associated  with  the  reference  coordinate  system. 
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BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


Specific  aims:  To  determine  the  single  and  double  scatter  cross  sections  for  random 
rough  surfaces. 

Results:  During  the  reporting  period  January  1  -  December  31,  1993,  the  principal 
investigator  submitted  four  papers  for  publication  in  Technical/Scientific  journals 
(see  item  #7).  One  has  been  accepted  for  publication,  two  others  have  been  revised 
and  resubmitted  and  the  fourth  is  still  under  review. 

In  the  paper  on  full  wave  solutions  for  rough  surface  bistatic  radar  cross  sec¬ 
tions,  item  #7(i)  extensive  comparisions  have  been  made  between  the  full  wave 
solutions  and  the  small  perturbation,  physical/geometric  optics  solutions  as  well  as 
published  numerical  and  experimental  results  for  horizontally  polarized  waves. 

In  the  second  paper,  item  #7(ii)  the  same  extensive  comparisons  are  made  for 
the  vertically  polarized  waves. 

In  the  third  paper,  item  #7(iii),  vertically  and  horizontally  polarized  diffuse 
double  scatter  cross  sections  have  been  evaluated  for  our  dimensional  random  rough 
surfaces  that  exhibit  enhanced  backscatter.  These  solutions  are  compared  with 
published  experimental  results.  A  comprehensive  parametric  study  is  also  conducted 
to  determine  the  dependence  of  the  level  and  the  angular  width  of  the  sharp  enhanced 
backscatter  upon  the  rough  surface  characteristics  (mean  square  height,  mean  square 
slope  and  the  mean  width  of  a  typical  depression).  The  results  for  the  double  scatter 
cross  sections  are  shown  to  depend  upon  the  polarization  of  the  electromagnetic 
waves.  This  is  consistant  with  the  experimental  results. 

The  last  paper,  item  #7(iv)  deals  with  enhanced  backscatter  from  random 
rough  surfaces  using  high  frequency-stationary  phase  approximations  to  the  full 
wave  solutions.  These  results  are  simpler  to  evaluate  and  provide  useful  physical 
insight  to  the  double  scatter  problem.  However,  the  angular  width  of  the  sharp 
enhanced  backscatter  cross  section  is  significantly  smaller  than  corresponding  results 
without  the  high  frequency  approximations.  Furthermore,  the  stationary  phase 
approximations  are  practically  independent  of  polarization. 

Plans  for  coming  year:  Extend  the  work  on  the  evaluation  of  the  single  and  double 
scatter  cross  section  from  one  dimenstional  random  rough  surfaces  to  two  dimen¬ 
sional  rough  surfaces. 
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